home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / dlantp.z / dlantp
Encoding:
Text File  |  2002-10-03  |  4.3 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAANNNNTTTTPPPP((((3333SSSS))))                                                          DDDDLLLLAAAANNNNTTTTPPPP((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLANTP - return the value of the one norm, or the Frobenius norm, or the
  10.      infinity norm, or the element of largest absolute value of a triangular
  11.      matrix A, supplied in packed form
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      DOUBLE PRECISION FUNCTION DLANTP( NORM, UPLO, DIAG, N, AP, WORK )
  15.  
  16.          CHARACTER    DIAG, NORM, UPLO
  17.  
  18.          INTEGER      N
  19.  
  20.          DOUBLE       PRECISION AP( * ), WORK( * )
  21.  
  22. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  23.      These routines are part of the SCSL Scientific Library and can be loaded
  24.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  25.      directs the linker to use the multi-processor version of the library.
  26.  
  27.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  28.      4 bytes (32 bits). Another version of SCSL is available in which integers
  29.      are 8 bytes (64 bits).  This version allows the user access to larger
  30.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  31.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  32.      only one of the two versions; 4-byte integer and 8-byte integer library
  33.      calls cannot be mixed.
  34.  
  35. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  36.      DLANTP returns the value of the one norm, or the Frobenius norm, or the
  37.      infinity norm, or the element of largest absolute value of a triangular
  38.      matrix A, supplied in packed form.
  39.  
  40. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  41.      DLANTP returns the value
  42.  
  43.         DLANTP = ( max(abs(A(i,j))), NORM = 'M' or 'm'
  44.                  (
  45.                  ( norm1(A),         NORM = '1', 'O' or 'o'
  46.                  (
  47.                  ( normI(A),         NORM = 'I' or 'i'
  48.                  (
  49.                  ( normF(A),         NORM = 'F', 'f', 'E' or 'e'
  50.  
  51.      where  norm1  denotes the  one norm of a matrix (maximum column sum),
  52.      normI  denotes the  infinity norm  of a matrix  (maximum row sum) and
  53.      normF  denotes the  Frobenius norm of a matrix (square root of sum of
  54.      squares).  Note that  max(abs(A(i,j)))  is not a  matrix norm.
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDLLLLAAAANNNNTTTTPPPP((((3333SSSS))))                                                          DDDDLLLLAAAANNNNTTTTPPPP((((3333SSSS))))
  71.  
  72.  
  73.  
  74. ARGUMENTS
  75.      NORM    (input) CHARACTER*1
  76.              Specifies the value to be returned in DLANTP as described above.
  77.  
  78.      UPLO    (input) CHARACTER*1
  79.              Specifies whether the matrix A is upper or lower triangular.  =
  80.              'U':  Upper triangular
  81.              = 'L':  Lower triangular
  82.  
  83.      DIAG    (input) CHARACTER*1
  84.              Specifies whether or not the matrix A is unit triangular.  = 'N':
  85.              Non-unit triangular
  86.              = 'U':  Unit triangular
  87.  
  88.      N       (input) INTEGER
  89.              The order of the matrix A.  N >= 0.  When N = 0, DLANTP is set to
  90.              zero.
  91.  
  92.      AP      (input) DOUBLE PRECISION array, dimension (N*(N+1)/2)
  93.              The upper or lower triangular matrix A, packed columnwise in a
  94.              linear array.  The j-th column of A is stored in the array AP as
  95.              follows:  if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
  96.              if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = A(i,j) for j<=i<=n.  Note
  97.              that when DIAG = 'U', the elements of the array AP corresponding
  98.              to the diagonal elements of the matrix A are not referenced, but
  99.              are assumed to be one.
  100.  
  101.      WORK    (workspace) DOUBLE PRECISION array, dimension (LWORK),
  102.              where LWORK >= N when NORM = 'I'; otherwise, WORK is not
  103.              referenced.
  104.  
  105. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  106.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  107.  
  108.      This man page is available only online.
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.